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A multiband empirical tight-binding model for group-III-nitride semiconductors with a wurtzite structure has 
been developed and applied to both bulk systems and embedded quantum dots. As a minimal basis set we 
assume one i-orbital and three p-orbitals, localized in the unit cell of the hexagonal Bravais lattice, from which 
one conduction band and three valence bands are formed. Non-vanishing matrix elements up to second nearest 
neighbors are taken into account. These matrix elements are determined so that the resulting tight-binding band 
structure reproduces the known F-point parameters, which are also used in recent k ■ p-treatments. Furthermore, 
the tight-binding band structure can also be fitted to the band energies at other special symmetry points of the 
Brillouin zone boundary, known from experiment or from first-principle calculations. In this paper, we describe 
details of the parametrization and present the resulting tight-binding band structures of bulk GaN, AIN, and 
InN with a wurtzite structure. As a first application to nanostructures, we present results for the single-particle 
electronic properties of lens-shaped InN quantum dots embedded in a GaN matrix. 
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I. INTRODUCTION 

Due to their unique physical properties, zero-dimensional 
semiconductor nanostructures, realized by either epitaxial 
growth or colloidal chemical synthesis,^ offer a broad range of 
applications.^ The three-dimensional confinement of spatially 
localized charge carriers in such tailor-made systems leads to 
a discrete and tunable one-particle spectrum, which can be 
used for a variety of optoelectronic applications, for quantum 
computing and quantum cryptography, and even for nanobio- 
logical applications like biological flourescence labeling.-^^ 

Most binary group II- VI and III- V semiconductor materials 
and their ternary and quaternary alloys crystallize in the cubic 
zincblende or in the hexagonal wurtzite phase, so the corre- 
sponding nanostructures can also be attributed to one of these 
two structures. Dependent on the material system and the ad- 
equate parameter range of the experimental conditions (e.g. 
growth temperature and substrate type for epitaxial growth, 
additionally the chemical environment and particle size for 
colloidal synthesis), it is even possible nowadays to realize 
either of the two structures for the same compounds. For ex- 
ample, epitaxially grown GaN/AlN quantum dots can be pro- 
duced in the metastable zincblende modification and in the 
thermodynamically stable wurtzite configuration.^ 

The calculation of the optical properties of such systems 
requires the knowledge of a set of single particle eigenstates 
and eigenvalues for the confined carriers (electrons and holes), 
which can be obtained by means of different methods. Rather 
simple models like effective mass approximations'^** can give 
a first insight into the behaviour of such systems. Multiband 
k • p-models^ri^ incorporate higher effects like valence band 
mixing, but still make use of the envelope function approxi- 
mation, thus not resolving the characteristic underlying lattice 
structure. Nevertheless, they have been succesfuUy applied to 
various material systems and were extended by the inclusion 
of strain and piezoelectricity effects. 

Empirical pseudopotential models (EPM)^^"'^ and empiri- 
cal tight-binding models (ETBM)'^"^^ allow for the possibil- 



ity of a microscopic description of nanostructures. While the 
EPM is capable of resolving variations on the atomic scale, it 
requires a large set of basis states, which limits the applica- 
tion of these models to small nanostructures. The ETBM uses 
a coarse graining on the scale of lattice sites, which makes 
it possible to stick to a small set of basis states and perform 
calculations on larger supercells with feasible effort. Addi- 
tionally, it gives a rather intuitive real-space picture of the 
system in terms of localized Wannier states. The coupling 
between different lattice sites is usually limited to first or sec- 
ond nearest neighbors, depending on the purpose. The goal 
is to analytically deduce a manageable set of equations for 
the tight-binding (TB) matrix elements in terms of bulk pa- 
rameters (e.g. the band gap, effective masses, spin-orbit split- 
ting) which can be accessed either from experiment or from 
first-principle calculations like DFT - LDA or recent GqWq- 
results.^^'^^ These TB matrix elements then enter the nanos- 
tructure calculation. 

When establishing an empirical tight-binding model, one 
can start from a Lowdin-orthogonalized atomic basis and use 
a linear combination of atomic orbitals (LCAO) as ansatz for 
the required eigenstates.^^ For nitride semiconductors with a 
wurtzite structure this has been done recently by Schulz et 
alP But it is as well justified to start from "effective bond- 
orbitals", i.e. Wannier-like orbitals localized within a unit 
cell. Within the LCAO spirit, these effective orbitals can in 
principle be expressed as linear combinations of the above 
mentioned atomic orbitals. Neither the atomic orbitals nor the 
effective orbitals are explicitly known or required within an 
empirical TB approach, as only the matrix elements between 
those orbitals are needed to obtain the TB band structure. 
Therefore, it is equally justified to perform the parametriza- 
tion directly for the effective bond-orbitals so that known bulk 
band structure properties are reproduced. Such a version of 
an ETBM is commonly called "effective bond-orbital model" 
(EBOM). The EBOM has the advantage that it usually al- 
lows for a better fit throughout the whole Brillouin zone (BZ) 
within a given basis set. 
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The EBOM has long been estabhshed for the cubic 
zincblende structure; a first EBOM parametrization by 
Chang^^ incorporated three-center overlap integrals in a basis 
set of one s- and three p-orbitals on each site of the fee Bra- 
vais lattice. This parametrization was restricted to coupling 
up to nearest neighbors, so that only the band energies at the 
F-point were fitted, besides the usual set of effective conduc- 
tion band masses and corresponding valence band parameters. 
Loehr augmented this model in Ref. by the inclusion of 
hopping up to second-nearest neighbors to additionaly fit the 
band-structure of the bulk material to an extended parameter 
set, including the X-point energies. This resulted in a better 
agreement of the resulting tight-binding conduction band with 
first-principle calculations.^^ 

To our knowledge, there exists only one parametrization 
of the EBOM for materials with wurtzite structure in the lit- 
erature.— As this work is restricted to a nearest neighbor 
parametrization and a fit to zone center energies only, we 
developed a new parametrization including second nearest 
neighbor matrix elements and a fit to band energies at other 
special BZ points. We apply this EBOM to the calculation 
of the electronic properties of lens-shaped InN quantum dots 
embedded within GaN. 

This work is organized as follows. In Sec. II, our specific 
EBOM is presented. We developed a second nearest neighbor 
parametrization. Furthermore, we describe the application to 
zero-dimensional nanostructures and discuss the inclusion of 
strain and piezoelectric fields. In Sec. Ill, the one-particle 
spectrum for a lens-shaped InN quantum dot embedded within 
GaN and a comparison with results from other k • p- and fully 
microscopic ETBM calculations is given. Section IV contains 
a summary, a conclusion and a brief outlook to possible ex- 
tensions of our model. 



II. THEORY 

A. Effective bond-orbital model for bulls semiconductors 

Linear combinations of atomic orbitals within one unit cell 
can be used as ansatz for the Wannier functions localized 
within the unit cell of the Bravais lattice. From these Wannier 
functions, the extended Bloch functions can be determined 
by means of a unitary transformation. As neither the atomic 
states nor the Wannier functions are expUcitly used in an em- 
pirical tight-binding model, it is not necessary to start from the 
atomic wave functions, but one can directly assume a basis of 
Wannier-like functions, which is the basic idea of the EBOM 
approach. 

As the conduction band wave functions at the BZ center 
predominantly transform i-like with some character, while 
the corresponding valence band wave functions transform like 
p-states with some s character, we use a localized sp^ basis 
per spin direction: 

|R,a), ae{s\,p,r\,pyt,Pzt,si,pU,Pyi,Pzi}- (1) 

Here R labels the sites of the hexagonal lattice, which is the 
underlying Bravais lattice of the wurtzite crystal structure. 



A trial wave function that satisfies the Bloch condition is 
the Bloch sum 

\yf^) = -Lycaik)l^e'^-^\R,a). (2) 

The band structure ^(k) is now given by the solution of the 
secular equation 

(k)c„.(k)=£(k)c„(k), (3) 

a' 

for each wave vector k, where 

//aa'(k)=IA(«-'*')£««;. (4) 
R,R' 

The EBOM matrix elements of the bulk Hamiltonian H^"^^ are 
thus given by 

= (R,a|//''""'|R',a'). (5) 

It should explicitly be pointed out that the artificial change of 
point group symmetry from C3,, (wurtzite) to Ce,. (hexagonal 
lattice) in the EBOM approach does not uniquely stem from 
the omission of the atomic basis, but rather from the specific 
set of basis functions used. For instance, the original inversion 
asymmetry of the wurtzite crystal could be restored when the 
set of basis functions, Eq. ([T]), is extended by states that are 
not parity eigenstates. This has been done for cubic systems 
by Cartoixa et al. in Ref. [s^- 

To include the influence of spin-orbit coupling, we follow 
Ref. [3TI As we expect the spin-orbit part of //bulk to be of 
weak influence, we assume only site-diagonal contributions, 
which stem from the p-orbitals. Additionally, the non-ideal 
c/a lattice constant ratio energetically seperates the p^- from 
the px- and the py-orbitals. These effects can properly be in- 
corporated by introduction of one spin-orbit splitting parame- 
ter Aso and one crystal field splitting parameter Aq-. 

When restricting the non-vanishing matrix elements, Eq. 
(|5]l, up to nearest or second nearest neighbors, the secular 
equation (|3]l can be solved analytically for high symmetry 
points throughout the BZ of the hexagonal lattice. This yields 
a set of equations for the EBOM matrix elements in terms of 
the energetic positions of the bands at the critical k-values. 
By expanding the elements of Eq. (|4]i around the BZ center 
and comparing the matrix representation to a corresponding 
k • p-Hamiltonian, ^^'^^ it is possible to deduce additional con- 
straints in terms of the conduction band effective masses and 
corresponding valence band parameters. 

The goal is to arrive at a solvable set of equations which link 
a sufficiently large number of to a desired set of band 
structure parameters. In practice, this will require the addi- 
tional omission of either matrix elements or band parameters, 
as the system of equations becomes rather complicated. The 
low symmetry of the hexagonal lattice will result in a larger 
number of independent parameters and equations than in the 
case of cubic crystal systems. An overview of the results for 
a coupling up to second nearest neighbors is given in Tab. U 
More details on the parametrization are given in App. [A] 
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TABLE I: Overview of the EBOM parametrization with coupling up 
to second nearest neighbors. The nomenclature for the band energies 
follows the usual single group notation, see e.g. Ref. 34. Note that 
an additional valence band parameter A-j has been neglected. 

Second nearest neighbor coupling: 26 EBOM matrix elements 



Band parameter 



Description 



,11 



me, 

Ai,A2,A3,A4,As,As 

Acr 

fC TV rv TV 

^1,3' -^1,3' ^2,4' '"1, 3' 
M\, Ml, Ml, M\ 

TjC TjV rjV 

, /I^ , lly 



direct band gap at F 
effective electron masses 
valence band parameters 
spin-orbit splitting 
crystal field splitting 
A-point energies 
L-point energies 
M-point energies 
/f-point energies 
Kane parameters 
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The band structure for this parametrization is now obtained 
by the diagonaUzation of the 8x8 matrix //„„/(k), Eq. (|4]l, 
for each k. 

For all calculations in the present paper, two distinct pa- 
rameter sets have been used. The first parameter set is de- 
rived from a consistent set of band parameters'^ obtained from 
Go Wo calculations based on exact-exchange optimized effec- 
tive potential ground states (OEPx),^^ supplemented by ad 
ditional band structure energies at high symmetry points 
Since the GoWo@OEPx band gaps and crystal field splittings 
still differ slightly from the experimental values, we decided 
to use a second parameter set, in which these values were re- 
placed by the parameters recommended by Vurgaftman and 
Meyer in 2003.^^ To obtain the correct band gaps, all conduc- 
tion band energies from the GqWo calculations were shifted by 
the respective difference in the second parameter set. More- 
over, we used the spin-orbit splittings of Ref. 37 in both sets, 
as the GqWo calculations did not include the electron spin. 
This should be a reasonable approach, as the spin-orbit split- 
ting is comparatively small in these systems. The two param- 
eter sets are listed in Tab. and will be referred to simply as 
"GqWo parameters" and "corrected GoWo parameters" from 
now on. 

In our opinion, the corrected GqWo parameters should 
clearly be preferred, as it is known that even highly sophisti- 
cated ab-initio approaches still do not properly reproduce the 
band gap. 

The resulting band structures are depicted in Fig. [T]for AIN, 
InN and GaN. The top of the valence band of each material 
is set to zero. One can easily identify the direct band gap 
in the Brillouin zone center, one spin-degenerate conduction 
band and three spin-degenerate valence bands, according to 
the employed basis set, Eq. ([T]i, of four orbitals per spin direc- 
tion. The twofold Kramers degeneracy of each energy level 
^(k) is a direct consequence of the time-reversal symmetry, 
as no external magnetic field is applied. Due to the fitting 
to the multiple high symmetry points on the BZ surface, each 



TABLE II: Empirical parameter sets used in the EBOM calculations. 
The first parameter set corresponds to results obtained in a DFT + 
GqWq treatment by Rinke et al.^'^-^^ while the second parameter set 
replaces some parameters by values recommended by Vurgaftman 
et at.?'' See the text for further discussion. Blank cells mean the 
adoption of the parameter of the alternate set. 
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-0.311 


-0.528 
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A3 


3.671 


5.414 


15.251 
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-1.147 
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As 


-1.329 


-2.510 


-7.060 










-1.952 


-3.202 


-10.078 
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8.844 
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3.355 


8.631 


5.972 


3.441 


AU [eV] 


-0.686 


-0.597 


-0.509 








A',' 3 [eV] 


-3.573 


-4.110 


-3.581 
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7.545 
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4.356 


7.332 
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4.442 


L'l 3 [eV] 


-1.515 


-2.065 


-1.732 








Ll ^ [eV] 


-1.689 


-2.144 


-1.838 








L\ y [eV] 


-6.033 


-6.984 


-5.769 








Ml [eV] 


8.084 


6.550 


4.934 


7.870 


6.821 


5.020 
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-0.837 


-1.111 


-0.997 








Ml [eV] 


-1.893 


-2.382 


-1.889 








M\ [eV] 


-3.649 


-4.518 


-3.714 








HI [eV] 


9.774 


7.982 


6.281 


9.560 


8.253 


6.367 


HI [eV] 


-0.914 


-1.609 


-I.40I 








HI, [eV] 


-5.202 


-6.474 


-5.422 









"For a given band gap, Ef,' are not independent parameters when are 
known (see Eqs. lAll - )A3t ) in App. |A] To obtain a better fit to the GoWq- 
band structure, these parameters can be adjusted by least-square fit values, 
see Ref. 24. When the band gap is subsequently altered, as in the set to the 
right, the analytic expression has to be used again. 



band has a finite bandwidth of a realistic magnitude, which the 
k • p-theory, of course, does not reproduce, as it is restricted to 
the vicinity of the BZ center within this basis set. In addition, 
no erroneous curvature of the bands into the band gap occurs 
for larger |k|. 

At first glance, the band structures do not differ signifi- 
cantly for both parameter sets. To emphasize the differences. 
Fig. 12] shows the band structure of InN around the F-point. 
By having a closer look, one can see that the energetic posi- 
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*^o^o parameters 

Corrected GqW^ parameters 




FIG. 2: (Color online) EBOM band structures for InN around T for 
the two parameter sets. Further details are again given in Tab. Ulland 
in caption of Fig. [T] 



tions are different, because of the different crystal field split- 
ting and band gap. Also, the curvatures of the bands differ for 
the two parameter sets. This is a result of the slightly different 
eI"^ and will be adressed again in Sec. HUB I More sophis- 
ticated methods for band structure calculation will give more 
conduction and valence bands in the energetic range around 
the band gap. Although this feature could also be included in 
our EBOM approach by augmenting the number of orbitals 
per unit cell, it would not only result in a more complicated 
parametrization, but also lead to significantly higher compu- 
tational costs for nanostructure calculations. Thus we stick to 
a minimal basis set of four bands per spin direction, which 
gives a reasonable agreement with the "true" band structure in 
the region of interest. The reliability of this basis set has also 



been established by the variety of existing eight-band-k • p cal- 
culations for these systems and comparisons to experimental 
results for device applications. ^^'^^ 

Of course, our EBOM parametrization is not limited to a 
specific set of parameters. The fit to the energies at the BZ 
boundaries allows for an adaption to a wide range of param- 
eters, as these additional constraints practically prevent spu- 
rious solutions where the bands curve into the band gap far 
away from the BZ center. These problems are widely known 
to occur for simpler k • p- and tight-binding parametrizations 
when an inappropriate set of F-point parameters is used. 

With our model, a more systematic investigation of the 
influence of single parameters on the properties of low- 
dimensional systems is possible, as the band structure is not 
sensitive to small perturbations in the input parameters. This 
stability of the parametrization transfers directly to the appli- 
cation on nanostructures, as e.g. spurios solutions in the bulk 
band gap will lead to coiTesponding states in the forbidden 
energy region of the nanostructure. In spite of the progress 
in the field of both sophisticated ab-initio calculations and 
highly refined experiments, this is and will remain an impor- 
tant feature, as certain physical quantities like band offsets, 
Luttinger parameters and optical matrix elements remain am- 
biguous because they are only indirectly measurable and de- 
pend on model assumptions. 



B. Application of the EBOM to quantum dots 

As we now have determined the EBOM matrix elements for 
the bulk materials, they can be used as input in the calculations 
for a quantum-confined nanostructure. In case of a quantum 
dot, the translational invariance is lost in all three spatial di- 
mensions, so the adequate ansatz for an eigenstate, Eq. (|2|i, 
is reduced to a direct linear combination of localized effective 
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orbitals 

|V/) = l^cR„|R,a). (6) 

a.R 

The corresponding secular equation is now given by 

H ^aa' CR'a' = Cro , (7) 
a',R' 

where E^^, are the EBOM matrix elements from Eq. ^ and 
the site indices R, R' now range over the finite sites of a suf- 
ficiently large supercell. In our present ip^-basis, the eigen- 
states and eigenenergies of Eq. O are obtained as the solu- 
tions of a 8A^ X 8A^ matrix eigenvalue problem. According to 
Refs.iliand[li a nanostructure made of one material A em- 
bedded in a barrier of material B can be modelled by using 
the matrix elements of the A-material for the corresponding 
lattice sites and vice versa. For the interface, a linear interpo- 
lation of the corresponding hopping matrix elements is used. 
The confinement potential for the carriers can properly be in- 
corporated by an upward shifting of the diagonal elements of 
the A-material by the valence band offset Afi,, between the 
two materials. When the band gap between the materials B 
and A exceeds this offset, we are naturally left with a type-I 
confinement potential for the electrons and holes. 

The application to one- or two-dimensional structures is a 
trivial task and can be done correspondingly. 

C. Possible inclusion of piezoelectricity and strain 

In the systems under consideration, there is always a spon- 
taneous polarization due to the deviation of the c/a-ratio from 
the value in the ideal wurtzite structure. Additionally, strain 
fields will influence the electronic properties of these sys- 
tems, if present. While there are in fact zero-dimensional sys- 
tems, like fully relaxed nanocrystals, where this effect can be 
neglected, the below presented model system of epitaxially 
grown InN quantum dots embedded in a GaN matrix will in 
fact be strained due to the lattice mismatch, and this strain will 
not only shift the band edges, but also alter the equilibrium 
positions of the lattice sites and thus the piezoelectric charge 
density. 

Both, the spontaneous and the strain induced polarization 
can be incorporated into the tight-binding calculations by the 
solution of the Poisson equation. '^'^'^ Although it has been 
discussed in previous publications like Ref. J_9, that for this 
specific quantum dot system a proper inclusion of a constant 
band-edge shift might be sufficient, a more general approach 
is of course desirable. 

Again, the one-to-one correspondence to the k • p model at 
r allows for a straight-forward inclusion of strain effects on 
the bulk band structure by augmenting the k • p-Hamiltonian 
by a strain-dependant part, as done in Refs. |32| and|33l The 
then obtained analytical dependance of the EBOM matrix el- 
ements on the deformation potentials a 12 of the con- 
duction band, Z)i_6 of the valence bands and the elastic stiff- 
ness constants can then be used either to determine a distance- 



dependant scaling law for the E^^,' (similar to the famous Har- 
rison d^^ ansatz^^) or directly be incorporated into the nanos- 
tructure Hamiltonian. In both cases, additionally an appropri- 
ate strain field has to be calculated for the low-dimensional 
system under consideration, either by atomistic^'^^ or contin- 
uum mechanical^^ approaches. 

As this extension of the EBOM requires a careful compari- 
son to further experimental and theoretical results, it is a topic 
of its own and part of ongoing research. Therefore, it will not 
furtherly be adressed in the present publication. In the fol- 
lowing section, we will neglect the influence of piezoelectric- 
ity and strain, in order to focus on the direct influence of the 
slightly different Kane parameters Ep''^ on the single particle 
results. 



ni. RESULTS FOR QUANTUM DOTS 

A. Model quantum dot geometry 

Earlier ETBM calculations of Refs. [HllOfor flie InN/GaN 
material system were performed using cubic supercells and 
fixed boundary conditions. In this paper we present a dif- 
ferent and improved kind of supercell, which is depicted in 
Fig. [3] by means of keeping the point group symmetry of 
the underlying hexagonal Bravais lattice in combination with 
periodic boundary conditions. This leads to several benefits, 
e.g. no artificial surface states can arise in the single-particle 
spectrum, and in contrast to cubic supercells our hexagonal 
one does not interfere with the Cg,, point group symmetry of 
the lattice. The simulated lens-shaped InN quantum dot has 
a diameter of 7.7 nm and a height of 3.1 nm. It is placed 
on top of a wetting-layer with thickness of one c lattice con- 
stant, since Stranski-Krastanov growth-mode is assumed for 
the given structure. The surrounding GaN supercell has a 
dimension of 36.4a x 42a x 14c with respect to the carte- 
sian axes. With this size, convergence for the one-particle 
wave functions is ensured. Furthermore, a completely strained 
structure is supposed, so that no deviations from the ideal lat- 
tice positions emerge. As valence band offset, we use the 
value recommended by Vurgaftman et alM- of AZi,, = 0.5 eV 
for both parameter sets. 



B. One-particle spectrum for embedded InN quantum dot 

The numerical diagonalization of the corresponding nanos- 
tructure Hamiltonian (using the folded spectrum method^'^) 
gives the desired single-particle states and eigenenergies 
around the energy gap of the quantum dot. We solve Eq. 
(|7]l for eight bound electron and hole states, using the EBOM 
parametrization for second nearest neighbors and taking spin- 
orbit coupling and crystal field splitting into account. 

The resulting eigenfunctions are visualized in Fig. |4]by 
isosurfaces of the probability density, supplemented by the 
respective eigenenergies. All states are invariant under rota- 
tions by ^ around the growth direction, according to the C^v 
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X-axis [nm] ya^is [nm] 



FIG. 3: (Color online) Geometry for the lens-shaped InN quantum 
dot on an InN wetting layer, embedded in GaN. The InN lattice sites 
of the QD are depicted with red dots, additionally, the intersections 
of the wetting-layer with the boundaries of the hexagonal supercell 
are visualized by the dotted red lines. 
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TABLE III: Orbital contributions of the sp^ basis to the bound one- 
particle states, calculated by summation over the site and the spin 
index of the projections (R, 0!|i/). Dominant parts are marked in 
bold. 



point group symmetry of the Bravais lattice. Each state is once 
again twofold degenerate due to time-reversal symmetry and 
well localized within the InN quantum dot. Tab. Hill reveals 
that the electron states mainly stem from the s-like conduction 
band, so a classification by their nodal structure is possible, 
y/j is 5-like, while i/l ™d Wi ^re complex linear combinations 
of the form p± = ^ [p^ ± ipy). y/| is a p^-like state, but dis- 
torted by the shape of the quantum dot. The hole states show 
similar transformation properties at first glance. Nevertheless, 
Tab. [Ill] reveals that at least two atomic p-states contribute 
to their formation, so that they underlie strong band-mixing 
effects. This is in agreement with results from several other 
multiband calculations-^^iffi 

By taking a closer look at the degeneracies in Fig. |4] we 
notice no fourfold degeneracy of i/l ^"d i/l, in contrast to 
k • p-calculations from Ref. |32l but reproduce the findings of 
earlier ETBM calculations from Ref. 2lL showing a twofold 
degeneracy for each bound state. The latter reference also in- 
cludes a detailed group theoretical dicussion of this issue. The 
only deviations we can report on is the fact that the EBOM 



eigenenergies of the electrons are more strongly bound for 
both parameter sets, compared to other ETBM results from 
Refs. [l^ and This discrepancy can safely be attributed 
to the once again different set of input parameters and thus 
does not contradict comparative studies from Refs. .41 and .421 
for zincblende nanostructures, which show a very good agree- 
ment between these approaches. Our results for the hole en- 
ergies agree well with the earlier calculations of Refs. [l9l a nd 
I20I if the piezoelectric field is neglected (see e.g. Ref. [l^ for 
details). 

A striking feature is the different order of the hole levels 
when switching between the parameter sets. The torus-shaped 
probability density belongs to the hole ground state t/zj' with 
the GqWo parameters, but is found as the first excited hole 
state 1/2 when using the coiTected GqWo parameters. Another 
look at Tab. |ll] reveals that the band gap of the dot material 
differs by less than 90 meV between these parameter sets; the 
confinement potential for the holes is even identical in both 
cases, as the same valence band offset is used. The crystal 
field splitting only differs by ca. 20 meV for both materials. 
Obviously, this rather small variation of the bulk band gap and 
the crystal field splitting, which also results in slightly differ- 
ent Kane parameters isj]'^ (see App. |A]l, suffices to change the 
level structure. Further studies (not shown) reveal that even 
variations in the order of magnitude of the accuracy of the 
input parameters can change the order. In addition, the differ- 
ing bulk band gaps of the two parameter sets lead to slightly 
different one-particle energy gaps Ef^ ~ e\—hi of ~ 0.97 
eV (Go Wo parameters) and ~ 1.05 eV (corrected Go Wo pa- 
rameters), respectively. When calculating optical properties 
like the excitonic absorption spectrum from the tight-binding 
single-particle spectrum, as e.g. done in Refs. Il9| and |42i 
the level structure and E^'^ are important characteristics. A 
change in the first will give rise to a change of the respective 
dipole matrix elements between the electron and hole states 
and thus alter the line intensity, while a change in Ef^ di- 
rectly shifts the energetic position of the line itself. 

As the variation of the dot size and the proper implementa- 
tion of strain effects and electrostatic built-in fields can addi- 
tionally alter the level ordering, the influence of different ma- 
terial parameters should be carefully investigated when dis- 
cussing the optical selection rules for a given geometry. For 
the system under consideration, we recommand the use of the 
corrected GqWo parameters. 



IV. CONCLUSION AND OUTLOOK 

In this paper, we have presented a multiband empirical 
tight-binding parametrization for both bulk semiconductors 
and nanostructures with a wurtzite structure, which represents 
an adaption of the effective bond-orbital model (EBOM) for 
the hexagonal phase. A basis set of one s- and three p-orbitals 
for each spin direction is placed on the sites of the underly- 
ing hexagonal Bravais lattice. Coupling up to second nearest 
neighbors has been used to fit one conduction and three va- 
lence bands to the energies and curvatures at the F-point and. 
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Go Wo parameters 



Vi ¥2 n>l n>l 

ei = 1 .4440 eV e2 = 1 -5802 eV 63 = 1 .5803 eV 64 = 1 .6858 eV 




n>2 n>3 n>4 



/!l= 0.4709 eV /!2 = 0.4676 eV /J3 = 0.4642 eV /J4 = 0.4544 eV 




Corrected Go Wo parameters 




FIG. 4: (Color online) Visualization of the probability density by isosurfaces of 75% (red), 45% (green) and 15% (blue) of the maximum value 
for electron \]/f and hole \j/f' states within the hexagonal supercell. The red dotted lines give the intersection of the wetting layer with the cell 
boundary. In addition, the corresponding one-particle eigenenergies are shown. All energies are given with respect to the valence band edge 
of GaN. The two parameter sets have been presented in Tab. |ll] Please note the interchanging of the states hi and /z2 when using the latter set 
of parameters. 



additionally, to the energies at high symmetry points through- 
out the whole Brillouin zone. The resulting band structures 
of the Ill-V-compounds InN, GaN and AIN were shown for 
two disctinct parameter sets, namely the newest GqWo- re- 
sults by Rinke et al. and a slightly modified set, in which 
we adjusted single critical parameters by replacing them by 
the values given by Vurgaftman et al. in order to obtain better 
agreement to experimental results. 

In addition, we demonstrated the application of this 
parametrization to low-dimensional structures. A lens-shaped 



InN quantum dot on an InN wetting layer, embedded in a GaN 
matrix, has been modelled within a hexagonally shaped super- 
cell with periodic boundary conditions. We have compared 
the resulting one-particle spectrum and the corresponding 
eigenstates to previous tight-binding results and have found 
a good concordance within the framework of the respective 
model. Furthermore, we have found that the two parame- 
ter sets yield a different order of hole states for the given 
dot diameter, although the corresponding bulk band structures 
barely differ at first glance. We strongly approve a careful 
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review of the set of material parameters used in such calcu- 
lations. For the present InN/GaN quantum dot system, we 
recommend the use of the corrected Go Wo parameter set. 

Besides the application to other quantum dot systems, 
like GaN in AIN, or different geometries, like coupled QDs 
or spherical nanocrystals, the present parametrization can 
easily be applied to one-dimensional (quantum wires) or 
two-dimensional (quantum wells and superlattices) structures. 
Moreover, the effects of strain and piezoelectric built-in fields 
can be incorporated on different levels of sophistication, as 
suggested in the second section of this publication. 
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Appendix A: EBOM parametrization for the hexagonal lattice 

The analytical dependance of the parameters 

^INx - (Al) 

of the band gap Eg, the spin-orbit and crystal field splittings 

Aso,Acr and the effective masses mjl'^ is given by the follow- 
ing two equations)^ 

p2^J^(m_^\ 3£g(A,„ + Eg) + A,,(2A,„ + 3Eg) 
" 2mo J 2Aso + 3Eg 

(A2) 



p2 [3£g(A,„ + Eg) + A,,(2A,„ + 3Eg)] 

^ 2mo J A„.A,„ + 3A,,.£g + 2A,„£g + 3£2 ' 

(A3) 

The EBOM parametrization scheme with coupling up to sec- 
ond nearest neighbors gives a set of equations which link the 
parameters of Tab. |II]to the EBOM matrix elements E^^, of 
Eq. |5] To obtain the desired number of free parameters for a 
one-to-one correspondance, one has to apply an adequate de- 
composition of the E^^, into two- and three-center integrals, 
following the guidelines of Ref. '25'. As the explicit solution is 
straightforward, but very unhandy in print, it shall not be given 
here in full form. Instead, further details will be made accessi- 
ble as supplementary material to this publication in mathemat- 
ical notation in Ref. 43 and, additionally, the explicit solution 
as MATLAB-compatible pseudocode in Ref. |4J, so that it can 
easily be used for own computations. 

In order to give at least a brief insight into the physical 
meaning of the E^, , we will give the results of the expan- 
sion of Eq. dUi to second order in k. In this limit, the EBOM 
and the k • p-presentation become equivalent. The following 
set of equations gives the EBOM matrix elements in terms of 
the parameters that were used in the 8-band-k • p-Hamiltonian 
of Refs. Il^andllll where the A, are Luttinger-like parameters 
which are connected to the anisotropic effective valence band 
masses. For the sake of simplicity, the parameter A7 has been 
set to zero in our approach. Its influence has turned out to be 
negligible.'*^ The upper index in now denotes R' — R 

in units of half the lattice constants a or c, respectively, so that 



, ka la mc 

R R = + y ey + — Cz. 



h ^1 

2ml 

h pI 

2mj- Eg 

h Pi 

2mj- Eg 

iP\\ 



_p(0,2,0) _ p(v^,l,2) _ ^p(0,2,2) _ 1 p(v/3,L0) _ . p(0,4,0) _ 9 (v/3,3.0) _ , (2v/3,2,0) \ 2 
_ 3^(^.1.0) _ 3^(V3.1.2) _ 3^(^3.0) _.^(2v/3,2.0) ^^{2^fifi)\ ^2 

2 2 ' ' / ' 



(2iEii 



2iE 



= {2iE, 



{-/iAfi) , 4.^jy3,l,2)^2/£(2v/3,0,0)^^.p(2v/3,2,0) 
(0,0,2) , g.^(y3,l,2) , ^..p (0,2,2) 



4/£L2^'2.0)+2/£i^-3-°)).«V3, 



(v/3.1,0) 



'2iE 



+ AiE, 



(0,2,2) 



-AiE 



(2^,2,0) ^g (v^,3,0) 



-AiE'^ 



(0.4.0) 



9 



A2+A4+A5- 



A2+A4-A5 
Ai+Aa 
2A, + 1 



A2+A4-A5 



A2+A4+A5- 



V2A, 



A1+A3 



A2 
A2 



Ai + ii 



_ ^1 



2 



,(^3,1.0) , (v/3.L2)_3 (^3,3.0) , (2x/3,2,0) _ „p(2v/3.0.0)A 2 

.(0,2,0) _ p(v^,l,2) _ ^p(0,2,2) _ 1 p(v/3,1.0) _ . p(0,4,0) _ 9 (v/3,3,0) _ , (2v/3,2,0)\ 9 
XX ^xx ^^xx 2 ^^xx 2 -^-^ -^^xx I : 

v^,l,2) _ ^ (0,2,2) (0,0,4) _ p(0,0,2)\ 2 

X -^^xx ^^xx ^xx / ' 

3 p(v/3,1.0) _ o p(x/5.1,2) _ 3 (^3,0) _ , (2x/3,2,0) _ , (2\/5,0,0)\ 2 

p(0,2,0) _ p(v^,l,2) _ ^p(0,2,2) _ 1 p(v/3,1.0) _ . p(0,4,0) _ 9 (v/3,3,0) _ ^ (2v/3,2,0)\ 2 

yy yy ^r^yy 2 >'>' 2 ^^yy I j 

4^1^3,1,2) _ 2£j0.2.2) _ 4£j0'''''*) _ £_(0,0,2)^ ^ ^2^ 

3^(v/3,1.0) _ 3£(y5.1.2) _ ^^1^3,3,0) _ g^p\/3,2,0) _ 2^(2x73,0,0)^ ^ ^2 

p(0,2,0) _ p(v^,l,2) _ ^p(0,2,2) _ 1 (\/3,l,0) _^p(0,4,0) _ 9 (^5.3.0) _ (2\/5,2,0)\ 2 
-'-'zz -'-'zz ^-'-'zz 2 zz 2 I J 

^^^(v^,i,2) _ 2£4'''^'^^ - 4£4°'''''*^ - £1°''''^'') • 



(A4) 



r 
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